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Quasi-power law ensembles are discussed from the perspective of nonex- 
tensive Tsallis distributions characterized by a nonextensive parameter q. 
A number of possible sources of such distributions are presented in more 
detail. It is further demonstrated that data suggest that nonextensive pa¬ 
rameters deduced from Tsallis distributions functions / (pt), Qi, and from 
multiplicity distributions (connected with Tsallis entropy), 92 , are not iden¬ 
tical and that they are connected via q\ + <72 = 2. It is also shown that 
Tsallis distributions can be obtained directly from Shannon information 
entropy, provided some special constraints are imposed. They are con¬ 
nected with the type of dynamical processes under consideration (additive 
or multiplicative). Finally, it is shown how a Tsallis distribution can accom¬ 
modate the log-oscillating behavior apparently seen in some multiparticle 
data. 
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1. Introduction 

The two most characteristic ensembles is, on one side, the one resulting 
in exponential (Boltzmann-Gibbs - (BG)) distributions, /g(X) ~ exp(— X/T), 
and, on the other side, the one with power distributions, fp(X) ~ A" -7 . 
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Both are encountered in the realm of the high energy multiparticle pro¬ 
duction processes investigated in hadronic and nuclear collisions. They are 
connected there with, respectively, nonperturbative soft dynamics operat¬ 
ing at small X’s and described by exponential distributions, /e(X)), and 
with perturbative hard dynamics, responsible for large X’s and described by 
power distributions, fp(X). These two types of dynamics are investigated 
separately. It is usually assumed that they operate in distinct parts of phase 
space of X , separated by X = Xq. However, it was found recently that the 
new high energy data covering the whole of phase space (cf., for example, 
mmm are best fitted by a simple, quasi-power law formula extrapolating 

between both ensembles mum-- 


H(X) = C ■ 




ex p(~Xl) for X—5-0, 
X~ n for X -»■ oo, 


(1) 


This formula coincides with the so called Tsallis nonextensive distribution 
[?] for n = l/(q — 1), 


h q (X) 


c n 


1-(1 -?)-£■ 
XqJ 


1-9 def n 

= C q exp q 




C i exp 



X 



This is the distribution we shall shall concentrate on and discuss. In Section 
[2] we shall discuss some examples of processes leading to such distributions. 
It depends on the nonextensivity parameter q and this can be different 
depending on whether it arises from a Tsallis distribution (qq) or from the 
nonextensive Tsallis entropy ( 52 )- Both are connected by qi+q 2 = 2 and this 
relation seems to be confirmed experimentally. This is presented in Section 
[H In Section [3] we shall discuss necessary conditions for obtaining a Tsallis 
distribution from Shannon information entropy. Section [5] demonstrates 
that a Tsallis distribution can also accommodate the log-periodic oscillations 
apparently observed in high energy data. Our conclusions and a summary 
are presented in Section [6l 


2. Some examples of mechanisms leading to Tsallis distributions 

In many practical applications, a Tsallis distribution is derived from 
Tsallis statistics based on his nonextensive entrop}Q, 


s q 


Pi ln g pi where ln g x 


pQ- 1 - 1 

q- 1 


( 3 ) 


1 d- IS] f° r most recent work with references; thermodynamical consistency of such 
an approach can be found in [Sj. We shall not discuss this point here. 
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On the other hand, there are even more numerous examples of physical 
situations not based on S q and still leading to quasi-power distributions in 
the Tsallis form. In what follows, we shall present some examples of such 
mechanisms, concentrating on those which allow for an interpretation of the 
parameter q. 


2.1. Superstatistics 

The first example is superstatistics [Ql (cf, also |ll| 12]) based on the 
property that a gamma-like fluctuation of the scale parameter in exponential 
distribution results in the g-exponential Tsallis distribution with q > 1 
(cf. Eq. Q). The parameter q defines the strength of such fluctuations, 
<7 = 1 + Var(X)/ < X > 2 . From the thermal perspective, it corresponds 
to a situation in which the heat bath is not homogeneous, but has different 
temperatures in different parts, which are fluctuating around some mean 
temperature Tq. It must be therefore described by two parameters: a mean 
temperature To and the mean strength of fluctuations, given by q. As shown 
in H3], this allows for further generalization to cases where one also has an 
energy transfer to/from heat bath. The scale T in the Tsallis distribution 
becomes then g-dependent: 

T = T eff = T 0 + (q-l)T v (4) 

Here the parameter Ty depends on the type of energy transfer, cf. mm 
for illustrative examples from, respectively, nuclear collisions and cosmic ray 
physics. 


2.2. Preferential attachment 

The second example is the preferential attachment approach (used in 
stochastic networks HSl). Here the system under consideration exhibits 
correlations of the preferential attachment type (like, for example, ”rich- 
get-richer” phenomenon in networks) and the scale parameter depends on 
the variable under consideration. If xo —> x' 0 (x) = xq + (q — l)x then the 
probability distribution function, /(x), is given by an equation the solution 
of which is a Tsallis distribution (again, with q > 1): 


df(x) 

dx 


*o( x ) 


/(*) 


/(*) = 


2 -q 
x 0 


1 - (1 - Q) 


x 


x 0 


1 

1-9 


( 5 ) 


For x' 0 (x) = xo one again gets the usual exponential distribution. 
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2.3. Multiplicative noise 

Consider now a Tsallis distribution from multiplicative noise mm- 
We start from the Langevin equation 


( 6 ) 


|+7 (t)p=m 


where q(t) and £(f) denote stochastic processes corresponding to, respec¬ 
tively, multiplicative and additive noises. This results in the following 
Fokker-Planck equation for the distribution function /, 


d / 

dt 


d{Kif) + d 2 (. K 2 f ) 


Stationary / satisfies 


with 


dp 

d{K 2 f) 

dp 


dp 2 


= KJ, 


(7) 

( 8 ) 


K\ = (£) — (j)p and K 2 = Var(£) — 2Cov{£ t ,^)p + Var(-f)p 2 . (9) 

In the case of no correlation between noises and no drift term due to additive 
noise (i.e., for Cov(£, 7 ) = (£) = 0 [IT]) its solution is a Tsallis distribution 
for p 2 , 

f(P ) = 


l + (9-l)f 


1 q , _ 2Var(£) 2V ardy) 

with T = - 7 ~zr ~~; -7 = 1 +- 77 -^. (10) 


<0 


(7) 


However, if we insist on a solution in the form of 


f(P ) = 


1 + 


P_ 

nT 


1 


n = 


Eq. (JHJ) has to be replaced by 

nT + p 


K 2 (p) = 


n 


Kiip) - 


q-T 

dl<\ jp) 
dp 


( 11 ) 


( 12 ) 


One then gets f{p) in the form of a Tsallis distribution, ( 1111 ) . but with 

(7) , , Var( 7 ) 


n = 2 + 


or q = 1 + 


(13) 


Var{ 7) (7) + 2Har(7) 

and with (/-dependent T (reminiscent of T e ff from Eq. (1141) discussed before 

cf. my- 


T{q) = (2 - q) [T 0 + (q - l)Ti] with T 0 = 


Covjf, 7) 

(7) : 


Ti = 


2(7)' 


(14) 
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2-4- All variables fixed 

Let us now remember that the usual situation in statistical physics is 
that out of three variables considered, energy U, multiplicity TV and tem¬ 
perature T, two are fixed and one fluctuates. Fluctuations are then given 
by gamma distributions m (in the case of multiplicity distributions where 
TV are integers, they become Poisson distributions) and only in the ther¬ 
modynamic limit (TV — > oo) does one get them in the form of Gaussian 
distributions, usually discussed in textbooks. In [19] we discussed in detail 
situations when two or all three variables fluctuate. If all are fixed we have 
a distribution of the type of 

f(E) = (1 - E/U) N ~ 2 3 . (15) 

This is nothing else but a Tsallis distribution with q = (TV — 3)/(TV — 2) < 

1 0 . 


2.5. Conditional probability 

For the constrained systems one gets q < 1. For example, if we have 
n independent energies, {-Ei=i,..., 7 v}) then each of them is distributed ac¬ 
cording to the Boltzman distribution, gt(Ei) = (1/A) exp (— Ei/X) (and 
their sum, E = F), is distributed according to gamma distribution, 

gN{E) = 1/[A(TV — l)](£'/A) Ar_1 exp(— E/X)). However, if the available en¬ 
ergy is limited, E = Not = const, then the resulting conditional probability 
becomes a Tsallis distribution with q < 1 0: 


f(Ei\E = Na) 


gi (Ei) g N _i (Na - Ei) (TV - 1) 


2 -q 


37V (TVa) 

Ei 


Na 


l-q 


TV — 3 

q= N^2 <1 ' 


A = 


aN 

TV- 1' 


]_E£\ 

N a ) 


N—2 


(16) 

(17) 


2 Actually, such distributions emerges directly from calculus of probability for situation 
known as induced partition [20]. In short: TV — 1 randomly chosen independent points 
{Ui,..., Un-i} breaks segment (0, U) into TV parts, length of which is distributed 
according to Eq. m- The length of the k th such part corresponds to the value of 
energy Ek = Uk+i — Uk (for ordered Uk). One could think of some analogy in physics 
to the case of random breaks of string in TV — 1 points in the energy space. Notice 
that induced partition differs from successive sampling from the uniform distribution, 
Ek € [ 0 , U — Ei — E2 — ... — Ek- 1], which results in f(E) = 1 /E [ 51 ] . 

3 One could get a Tsallis-like distribution with q > 1 only if the scale parameter A 
would fluctuate in the same way as in in the case of superstatistics, see [22] , 
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2.6. Statistical physics 

We end this part by a reminder of how Tsallis distribution with q < 1 
arises from statistical physics considerations. Consider an isolated system 
with energy U = const and with v degrees of freedom (n particles). Choose 
single degree of freedom with energy E (i.e., the remaining, or reservoir, 
energy is E r = U — E). If this degree of freedom is in a single, well de¬ 
fined, state then the number of states of the whole system is £l(U — E) 
and probability that the energy of the chosen degree of freedom is E is 
P(E) oc Q(U — E). Expanding (slowly varying) In 17(E) around U, 


In 17(1/-E) 


y, 1 d {k "> In 17 

k =0 


with /3 


1 def d In 17 (E r ) 
kgT dE r 


and (because E « U) keeping only the two first terms one gets 
lnE(E) oc In 17(E) oc — /3E, or P{E) oc exp(— /3E), 


(18) 


(19) 


i.e., a Boltzmann distribution (or q = 1). On the other hand, because one 
usually expects that 17 (E r ) oc {E r f v) OLlv ~° 12 (where aq^ are of the order of 
unity and we put a i = 1 and, to account for diminishing the number of 
states in the reservoir by one, «2 = 2) [23], one can write 


d k (3 

m k 


oc (-1 ) k k\ 


u-2 

pk+l 

H/f 


(-1 ) k k\ 


p k ~ l 

{y - 2) k 


( 20 ) 


and write the full series for probability of choosing energy E: 


P{E) oc ^Zrr E = ex P 


Q(U) 

= c( 1 - 


E 


(-E 


i 


i 


-BE 


k=0 k + l(n-2Y 
iy- 2 ) 


{-PE) 


fc+i 


v — 2 


( 21 ) 


(where we have used the equality ln(l+x) = —l) fc [a: fc_, “ 1 /(/r+1)]). This 

result, with q = 1 — 1 j[y — 2) < 1, coincides with results from conditional 
probability and the induced partition. 


2. 7. Fluctuations of multiplicity N 

Constant values of U. N and T result in Eq. m and q < 1. To get 
larger values of q , one has to allow for fluctuations of one of the variables: 
U. N or T. In superstatistics ns mug it was T that was fluctuaing, let 
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us now consider the example of fluctuating N. It means that, whereas for 
fixed N (to simplify notation we changed N — 2 to N here) 

/ E\N 

f N (E) = (1 — — J and U = E = const , (22) 

for N fluctuating according to some P(N ) the resulting distribution is 

f(E) = Y / fN(E)P(N). (23) 


The most characteristic for our purposes are situations provided by the, 
respectively, Binomial Distribution (BD), Ea. (l24D . Poissonian Distribution 
(PD), Ea. (l25l) and by the Negative Binomial Distributions (NBD), Ea. (l2fll) 

(cf. my 


Pbd(N) 
Ppd(N) 
Pnbd{N ) 


M\ ( <N> \ N (\ < n >\ M ~ N 

N\(M - N)\ V M ) \ M ) 


Nl 

T(N + k ) /<iV>\ iV / < N >\~ k ~ N 

T(N + l)r(fe) \ k J V 1+ k ) 


(24) 

(25) 

(26) 


They lead, respectively, to Tsallis distribution with q = 1 — 1/M < 1, 
Ea. (1271) . to exponential Boltzmann distribution, Eg. (1281) with q = 1, and to 
Tsallis distribution with q = 1 + 1/k > 1, Ea. (l29l) (/? = ( N)/U ): 


fp D (E ) = exp (~/3E), 

/™w-(i + f)-*. 

Note that in all three cases 

_ Var(N) 1 
Q ~ ~ < N > 2 ~ < N > 


(27) 

(28) 
(29) 


(30) 


It is natural that for BD where Var{N)/ < N >< 1 one has q < 1, for PD 
with Var(N)/ < N >= 1 also q = 1 and for NBD where Var(N)/ < N » 
1 one has q > 1 (cf., also case U = const and T = const considered in [19]). 

Fluctuations of N can be translated into fluctuations of T. Notice first 
that US] (cf. also [ 25 ]) NBD, P(N) = T(N + k)/[T(N + l)r(As)] • 7 fc (l + 
7 arises also if in the Poisson multiplicity distribution, P(N) = 
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N n e N /N\, one fluctuates the mean multiplicity N using gamma distribu¬ 
tion f(N) = r y k N k ~ 1 /T(k) ■ e -7 ^ with 7 = k/ < N > Q. Now, identifying 
fluctuations of mean N with fluctuations of T, one can express the above 
observation via fluctuations of temperature. Noticing that P = N/U (i.e., 
that < N >= U < j3 > and 7 = k/[U < P >]) one can rewrite f(N) as 


fiP) 


k ( k/3 \ k 1 f kp \ 

<p>m exp ( v _ 7j7j“ 


Vg-i <p>) 


(«-i)<^>r(;i I ) 



£-V 

q-l<P>) 


(31) 


And this is just a gamma distribution describing fluctuations of ft = 1/T 
discussed in HU. 

There is more to the physical meaning of q. Because U =< N > T 
the heat capacity C can be written as 1/C = dU/dT =< N >. However, 
because in our case U = const , i.e., Var(N)/(N) 2 = Var(f3)/(/3) 2 (or 
< N >~< P > and Var(N) ~ Var(p)), the formula q — 1 = Var(N)/ < 
N > 2 —1/ < N > obtained before coincides with a similar formula obtained 


in [25] : 


Var(p) l Var{T) 1 
<p > 2 ~C~ <T > 2 “ C‘ 


(32) 


This can be confronted with q — 1 = Var(N)/ < N > 2 from [13] . which, 
because of the form of f{fi) in Eq. (I3T1) . can be rewritten approximately as 


, Var(p) VariT) 

q 1 “ <I>^ ~ 7^' 


(33) 


3. Tsallis distribution from Shannon entropy 

Tsallis distribution is usually derived either from Tsallis entropy (via the 
MaxEnt variational approach, not discussed here) or from some dynamical 
considerations, some examples of which are presented in this paper. How¬ 
ever, it turns out that it also emerges in a quite natural (in the same Max¬ 
Ent approach) way from Shannon entropy, provided one imposes the right 
constraints. In fact, as shown in [28], one can establish a transformation 
between these two variational problems and shown that they contain the 

4 We have two type of averages here: X means average value in a given event whereas 
< X > denotes averages over events (or ensembles). 
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Fig. 1. (Color online) Examples of q < 1 or q > 1. (a) Distributions of 

p r T el = \pxp? et \/ \pj e t | for particles inside the jets with different values of Pj e t- 
Distributions are fitted using a Tsallis distribution ([2]) with T = 0.18. (b) 

(N + 1 )P{N + 1 )/P(N) = a + bN as function of multiplicity N in jets with dif¬ 
ferent values of pj e t as presented in (a). Depending on the range of phase space 
covered, the corresponding values of q are 0.885 for pj et G (4 — 6) Gev and 1.094 
for pj et G (24 — 40) GeV. Data are from [31 EH and from the Durham HepData 
Project; http://hepdata.cedar.ac.uk/view/irn9136932. 


same information. This means that the two approaches seem to be equiv¬ 
alent, one can either use Tsallis entropy with relatively simple constraints, 
or the Shannon entropy with rather complicated ones (cf., for example, a 
list of possible distributions one can get in this way [29]). 

In general, Shannon entropy for some probability density f(x), S = 
— f dxf(x)ln[f(x)], supplied with constraint < h(x) >= f dxf(x)h(x) = 
const, where h(x) is some function of x, subjected to the usual MaxEnt 
variational procedure, results in the following form of f(x): 

fix) = exp [A 0 + \h(x)] , (34) 


with constants Ao and A calculated from the normalization of f(x) and from 
the constraint equation. It is now straightforward to check that 

< z >= zo = - - where z = In 

2 -q 


, , E 

1-(1 - Q ) t F 
1 o. 


(35) 


results in f(z) = (1/zo) exp (— z/zq) which translates to (remembering that 
f{z)dz = f(E)dE) a Tsallis distribution 


m 


i 



Zp E 
1 + zq To 


1 + Z Q 

z 0 


2 ~q 


, E 

l-(l-g)- 

4 0 j 


1 

1 ~q 


■ (36) 


(l+^o) Tq 
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The parameter To can be deduced from the additional condition which must 
be imposed, namely from the assumed knowledge of the (. E ) (notice that in 
the case of BG distribution this would be the only condition). 

So far the physical significance of the constraint (13511 is not fully under¬ 
stood. Its form can be deduced from the idea of varying scale parameter in 
the form of the preferential attachment, Eq. ([5]). which in present notation 
means T —> T(E ) = To + (q — 1 )E. As shown in (J3J) it results in Tsallis 
distribution (13U1) . This suggest the use of z = In [T(E)/Tq\ constrained as in 
Eq. (1351) . In such approach In f(E) = — [1 /{q — l)]ln[T(E)] + [(2 — q)/(q — 
l)]ln(To) and, because S = —(In/(E)), therefore S = 1/(2 — q) -fin (To) 
for Tsallis distribution becoming S = 1 + ln (To) for Boltzmann-Gibbs (BG) 
distribution (q = 1). 

It is interesting that the constraint (13511 seems to be natural a for mul¬ 
tiplicative noise described by the Langevine equation: dp/dt + 'y(t)p = £(i), 
with traditional multiplicative noise y(t) and additive noise (stochastic pro¬ 
cesses) £(t)) (see [18] for details). In fact, there is a connection between 
the kind of noise in this process and the condition imposed in the MaxEnt 
approach. For processes described by an additive noise, dx/dt = £(t), the 
natural condition is that imposed on the arithmetic mean, < x >= c+ (£)t, 
and it results in the exponential distributions. For the multiplicative noise, 
dx/dt = x'y(t), the natural condition is that imposed on the geometric mean, 
< lnx >= c+ ( 7 )t, which results in a power law distribution [30]. 

4. Tsallis entropy vs Tsallis distributions (qi -f <72 = 2) 

One has to start with some explanatory remarks. The Tsallis distribu¬ 
tion can be also obtained via MaxEnt procedure from Tsallis entropy, 

Sq = = ( ln 9 Pi) q = -(^2-qPi)q=l, (37) 

where \n q x = ( x 1 ~ q — 1) /(I — q) and (x) q = Y/Pi x i- Now, depending on 
the condition imposed one gets from S q 

either f(x) = q[ 1 + (1 — q)x]~ T ~Q for (x)\ (38) 

or f(x) = (2 - q)[l + (q - l)x]^ for (x) q . (39) 

However, after replacement of q by q\ = 2 — q, the distribution (]38l) becomes 
the usual Tsallis distribution (|39D . Therefore one encounters an apparent 
puzzle, namely the q\ of Tsallis distribution does not coincides with the 
q 2 of corresponding Tsallis entropy, instead they are connected by relation 
qi + q 2 = 2. The natural question therefore arises: is such a relation seen in 
data? As shown in m that seems really be the case, at least quantitatively. 
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Fig. 2. (Color online) (a) Energy dependencies of the parameters q obtained from, 
respectively: multiplicity distributions P(N) [35] (squares), from different analysis 
of transverse momenta distributions / (pt) in P+P data ( [31] - circles, full symbols) 
and from data on / ( pp ) from Pb+Pb collisions ([25] - half filled circles). ( b ) Energy 
dependence of the charged multiplicity for nucleus-nucleus collisions divided by 
the superposition of multiplicities from proton-proton collisions fitted to data on 
multiplicity taken from [35] (NA49) and from compilation ,3Bj. 


This is seen when comparing q = q\ obtained from data on pp distributions 
(cf. Fig. [2(a)) to g = (72 obtained from data on multiplicities in p — A 
collisions assuming that entropy is proportional to the number of particles 
produced (cf. Fig. [2] (£>))- Whereas q\ is deduced from a Tsallis distribution 
taken in one of the forms discussed above, q 2 is deduced directly from the 
corresponding entropy S q of the p — A collision. Assume that such collision 
can be adequately described by a superposition model in which the main 
ingredients are v nucleons which have interacted at least once [32]. Assume 
further that they are identical and independent and produce ni secondaries 
of each other. As a result a N = Ya =i n i are produced in one collision and 
the mean multiplicity is (TV) = ( n)(rii ) (where (v) is the mean number of 
nucleons participating in the collision and ( ni) the mean multiplicity in i th 
elementary collision. The corresponding entropy Sq ^ of such process will 
then be g-sum of v entropies Sq 1 ^ of individual collisions and is given by: 


•S'"’ = £ 


v\ 


^ (y - k)\k\ 


(1-9) 


k -1 


S 


k 


1 + (1 - q)S± 


(1) 


-1 


1-9 


(40) 


Notice that In 


1 + (1 - q)Sq 




= z/ln 1 + (1 — q)Sq 1 ' 1 and Sq^ v■ . 

For q < 1 entropy S^ is nonextensive because S^ jv qq. For q > 1 
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one has S^ > 0 only for q < 1 + 1 / Sq 1 ' 1 and Sq ^/ v 0, i.e., entropy is 

extensive, 0 < Sq^/u < S*q\ 

Because ( Naa) > Np{N pp ) = v(N pp ), the nonextensivity parameter 
obtained from the corresponding entropies must be smaller than unity, q 2 < 
1. On the other hand, all estimates of nonextensivity parameter from Tsallis 
distributions lead to q\ > 1. 

5. Dressed Tsallis distributions (log-periodic oscillations) 

The pure power-like distributions are known to be in many cases deco¬ 
rated by specific log-periodic oscillations (i.e, multiplied by some dressing 
factor R ) [37]. They suggest some hierarchical fine-structure existing in the 
system under consideration and are usually regarded as possibly indicating 
some kind of multifractality in the system. Closer inspection of recent data 
from LHC [T11311211351 39] reveals that, for large transverse momenta pp. one 
observes a similar effect, cf. Fig. [3] So far, the prevailing opinion is that 
this is just an apparatus induced artifact with no meaning. However, its 
persistence in the type of experiment considered, energy or type of collision 
process (provided that the range of pp covered is large enough) calls for 
some explanation. Such an explanation was considered in [401141] . In [37] . 
the only possibility investigated was to attribute such oscillations to com¬ 
plex values of the power index n in Eq. m ■ As shown in m and also below, 
it also works in the case of quasi-power like Tsallis distributions. However, 
because one now also has a scale parameter T and a constant term, one can 
also offer another explanation: real n but log-periodically oscillating scale 
parameter T. This was discussed in detail in EH- We shall present both 
possibilities here. 


5.1. Complex nonextensivity parameter q 

For simple power laws, one has some function 0(x) which is scale invari¬ 
ant, O(Xx) = pO(x) and 0(x) = Cx~ m with m = — In pi In A. However, 
this can be written as pX m = 1 = e* 2?rfc , where k is an arbitrary integer. We 
then have not a single power m but rather a whole family of complex pow¬ 
ers, mfc, with mfc = — In In A + *27rfc/ln A. Their imaginary part signals 
a hierarchy of scales leading to log-periodic oscillations. This means that, 
in fact, 0(x) = J2k =o WkRe ( x~ mk ) = x~ Re( ^ mk ^ J2k =o w k cos V m {mk) hi(ic)] 
(where Wk are coefficients of the expansion). This the origin of the usual 
dressing factor appearing in m and used to describe data: 

R(E) = a + b cos [c\n(E + d) + f] (41) 

(only w\ and u >2 terms are kept). 
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Fig. 3. (Color online) Examples of log-periodic oscillations, (a) dN/dpr for the 
highest energy 7 TeV, the Tsallis behavior is evident. Only CMS data are shown 
PQ, others behave essentially in the identical manner. ( b ) Log-periodic oscillations 
showing up in different experimental data like Qj or ATLAS [3] taken at 7 TeV. 
(c) Results from CMS [T] for different energies, (d) Results for different systems 
(p + p collisions compared with Pb + Pb taken for 5 % centrality [38] . Results from 
ALICE are very similar [35]. 


It turns out that a similar scaling solution can also be obtained in case 
of a Tsallis quasi-power like distribution. To this end one must start from 
stochastic network approach, Section 12.21 and Eq. ([5]) , in which Tsallis 
distribution is obtained by introducing a scale parameter depending on the 
variable considered. In our case it is df(E)/dE = —f(E)/T(E) resulting in 


m 


n- If E \ ~ n 
nT 0 \ + nT 0 ) 


for T(E) 



n. 


(42) 
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In final difference form (with change in notation: Tq replaced by T) 


df(E) 

dE 


HE) 

T(E) 


f(E + 5E) 


—n5E + nT + E 
nT + E 


f(E). 


(43) 


We consider a situation in which 5E = anT(E) = a(nT + E). It depends 
now on the new scale parameter a (a < 1/n in order to keep changes of 5E 
to be of the order of T) and can be very small but always remains finite. 
It can now be shown that f[E + a(nT + E)] = (1 — an)f(E) which, when 
expressed in the new variable x = 1 + E/(nT), corresponds formally to the 
following scale invariant relation: 


g[( 1 + a)x\ = (1 - an)g(x). (44) 

Following the same procedure used to obtain dressed solutions of scale in¬ 
variant functions discussed at the beginning of this Section, one arrives at 
the dressed Tsallis distribution (we keep, as before, only the two lowest 
terms, k = 0 and k = ljf|: 


9(E) 


E \~ mo 

1 + — I u; 0 + wi cos 


2vr / E_ 

ln(l + a) n V nT 


(45) 


with mo = — ln(l — cm)/ ln( 1 + a) n. In addition to the scale parameter 
a one has two more parameters occurring in the dressing factor R, wq and 
w\. The other parameters occurring in Eq. ed are expressed by the original 
parameters in the following way: a/b = wq/w\, c = 27r/ln(l + a), d = nT 
and / = —27rln(nT)/ln(l + a) = —clnd. One can, however, consider a 
more involved evolution process, with k sequential cascades; in this case the 
additional parameter n changes parameter c in (ED, c —> d = c/k. It does 
not affect the slope parameter mo but changes the frequency of oscillations 
which now decrease as 1 /k. Comparison with data requires n 22 (cf. mi 
for details). 


5.2. Log-periodically oscillating T 

As mentioned before, one can translate a dressed Tsallis distribution 
into a normal one but with a log-periodically oscillating in pr scale factor 
T, cf. Fig. m. The formula used there to fit the obtained results resembles 
that for dressing factor ED, 

T = a + b sin [c (ln(Fl + d) + f] . (46) 

5 Notice that in Eq. m n ^ mo- However, n and T are both unknown a priori 
parameters. Therefore, for fitting purposes, where we have to use two (and not three) 
parameter Tsallis distribution, we use [1 + E/(moT')]~ rn ° with fitting parameters mo 
and T' . In terms of T, n and a we have T' = nT/mo = —nT ln(l + a)/ln(l — an). 
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Fig. 4. (Color online) (a) Oscillations of scale parameter T leading to identical 
dressed Tsallis distribution as shown in Fig. [3 Jd (obtained for CMS data at 7 TeV 
and fitted using Eq. lHSl) h ( b ) Dependencies of t/tq from Eq. (1541) and £ — £o from 
Eq. (I5TD resulting in oscillations of T shown in panel (a). 


In fit shown in Fig. |4h parameters (generally energy dependent) are a = 
0.143, b = 0.0045, c = 2.0, d = 2.0, / = -0.4. 

To explain Eq. (|46l) one uses a stochastic equation for the temperature 
evolution [42] written in Langevin formulation with energy dependent noise, 
£(t,E), and allowing for time dependent E = Eft) 


^ + 7 r + «‘' E)T = 4 - (47 > 

Assuming now a scenario of preferential attachment (cf. Section 12.21 above) 
known from the growth of networks IE]) one has 


dE E m 

-TT ~ - \-T. 

dt n 


and Eq. (1471) has now the form: 


(48) 


(- + t) % + -T + E)T = 4>. (49) 

V n J dE t 


After straightforward manipulations (cf. [40] for details) one gets, for large 
E (i.e., neglecting terms oc 1/E): 


1 d 2 T 
n d{\n E) 2 


+at,E) 


dT 


■d£(t,E ) 

+ T = 0. 


ci(ln E) d(ln E) 


(50) 


Notice the change of notation, we discuss formulas for energy E but results are for 
transverse momenta pr here. However, they are taken at the midrapidity, i.e., for 
y — 0, and for large transverse momenta, pr > M, and in this region one has 
E = yjM 2 + cosh(y) 
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Assuming now that noise £(f, E ) increases logarithmically with energy, 

,2 


UJ 


Z{t,E) = &(*) + — InE. 


n 


(51) 


In this case Eq. (15011 becomes an equation for the damped hadronic oscillator 
with solution in the form of log-periodic oscillation of temperature with 
frequency uj and depending on initial conditions phase shift parameter </>: 


T = C exp < — n 


1 gt,E) 

2 T 2 


In E > • sin(u; In E + (j)). 


(52) 


Averaging the noise fluctuations over time t and taking into account that the 
noise term cannot on average change the temperature, 1/r + (£( t,E )) = 0, 
one arrives at 


T = a + — sin(w In E + </>). 
n 


(53) 


This should now be compared with the parametrization of T(E) given by 
Eq. (1151) and used to fit data in FigH] 0. 

We close with the remark that, instead of using energy dependent noise 
£(t, E) given by Eq. (1511) and keeping the relaxation time r constant, we 
could equivalently keep the energy independent white noise, £(t,E) = £o (t), 
but allow for the energy dependent relaxation time, for example in the form 
of 


T = t(E) 


nT 0 

n T (jj^ In E 


(54) 


In this case the temperature evolution has the form 


T(t) = (T) + [T(t = 0) - (T)]E- tuj2 / n exp , (55) 

and T gradually approaches its equilibrium value (T). Actually, for r = 
t(E ), as in our case, this approach towards equilibrium is faster for large E. 
This is because, in addition to the usual exponential relaxation characteristic 
for r = const case, we have an additional factor ~ E~ tw ! n . 


6 . Summary and conclusions 

We presented examples of possible mechanisms resulting in quasi-power 
distributions exemplified by Tsallis distribution, Eq. (J5J). Our presentation 
had to be limited, therefore we did not touch thermodynamic connections 

7 Notice that only a small amount of T, of the order of b/a ~ 3%, emerges from the 
stochastic process with energy dependent noise; the main contribution comes from 
the usual energy-independent Gaussian white noise. 
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of this distribution UM or the possible connection of Tsallis distributions 
with QCD calculations discussed recently (33, 33]. 

The main results presented here can be summarized in the following 
points: 


• Statistical physics consideration, as well as ’’induced partition pro¬ 
cess”, results in Eq. (USD, i-e-, in Tsallis distribution with q = (IV — 
3)/(N — 2) < 1. Fluctuations of the multiplicity N modify the pa¬ 
rameter q which is now equal to q = 1 + Var(N)/(N) 2 — 1/(1V), cf. 
Eq. (1301) . Notice that conditional probability for the BG distribution 
again results in Eq. (fl5l) . 


• Fluctuations of the multiplicity N are equivalent to results of an ap¬ 
plication of superstatistics, where the convolution 

f(E) = J < 7 (T)exp(—(56) 

becomes a Tsallis distribution, Eq. (J2D , for 


g(T) 


1 fnT 0 \ n 

exp 



(57) 


• Differentiating Eq. (1561) one gets 

= where r(£) = r « + !• (58) 

This is nothing else than a ’’preferential attachment” case, again re¬ 
sulting in a Tsallis distribution which for T(E ) = Tq becomes a BG 
distribution, cf. Eq. ([5]). 

• Replacing in Eq. (1581) differentials by finite differences, cf. Eq. (1331) . 
one gets for 5E = anT(E) the scale invariant relation, Eq. (|44p . which 
results in log-periodic oscillations in Tsallis distributions @. 

In addition to this line of reasoning, we have also brought in the prob¬ 
lem of the apparent duality between the nonextensive parameters obtained 
from the whole phase space measurements of multiplicity and more local 

8 Among numerous other expianations we can therefore say that we have demonstrated 
that a Tsailis distribution, which can be regarded as generalization to real power n 
of such well known distributions as the Snedecor distribution (with n = (y + 2)/2 
with integer v, for v —\ oa it becomes an exponential distribution), can be extended 
to complex nonextensivity parameter. 
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measurements of transverse momenta. This point deserves an experimental 
and phenomenological scrutiny. 

Finally, we tentatively suggested that, by choosing the right constraints, 
which account for additive or multiplicative processes considered, one can 
also get a Tsallis distribution directly from the Shannon information en¬ 
tropy. 
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